Conducting Jahn- Teller Domain Walls in Undoped Manganites. 
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We investigate the electronic properties of multi-domain configurations in models for undoped 
manganites by means of variational and Monte Carlo techniques. These materials display simul- 
taneous Jahn- Teller distortions and magnetic ordering. We find that a band of electronic states 
appears associated with Jahn- Teller domain walls, and this band is localized in the direction per- 
pendicular to the walls. The energy and width of this band depends on the conformational properties 
of the domain walls. At finite temperatures, the conductance along the domain walls, induced by 
the localized domain wall bands, is orders magnitude larger than in the bulk. 
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I. INTRODUCTION 

In the quest to understand and control the proper- 
ties of broken symmetry phases, domain walls play a 
central role. Among the many materials with ground 
state broken symmetries, strongly correlated transition 
metal oxides, and manganites among them, are particu- 
larly interesting. In these compounds several degrees of 
freedom are simultaneously important, and different bro- 
ken symmetry phases with similar characteristic energy 
scales either compete, as in phase separated materialsji 
or coexist, as it occurs in multiferroicsi^ Research on 
magnetic domain walls in general has been particularly 
intense, and it has proven to be extremely important to 
explain the static and dynamic properties of magnetic 
materials. Early work^ led to the development of a va- 
riety of important concepts to understand domain walls. 
Technological applications, and the need to control and 
understand the important details of domain walls, stim- 
ulated considerable and wide research on this topic r'"— 
unveiling a variety of interesting aspects of these walls. 
In particular, the magnetic domain walls of manganites, 
the materials with the colossal magnetoresistance, have 
also attracted much attention. Already in Ref. 10, the 
large magnetoresistance was attributed to domain wall 
scattering, an hypothesis that led to both experimental 
and theoretical work on the subject of resistance of do- 
main wallsJ^"— Recently, it was shown that magnetic 
domain walls in a ferromagnetic metallic material could 
be insulating jiS. 

Similarly, domain walls and gradients of the order pa- 
rameter play a crucial role in our understanding of other 
collective phenomena, such as superconductivityj^*' and 
ferroelectrityj ^^i^^ Early theoretical work showed that 
physical properties absent in bulk materials can arise in 
domain walls, ^^ and the conductivity of ferroelastic do- 
main walls has been showed to be different from the con- 
ductivity in the bulki^^'-- Another very interesting field 
is the study of properties of domain walls in multifer- 
roics, or magnetoelectric materials. Different orderings 
can change across domain walls in these materials, and 



indeed it has been observed than often the change in 
one order parameter is correlated with modifications in 
another *2S, The interplay between order parameters can 
affect the physical properties of domain walls. In the 
antiferromagnetic and ferroelectric phases of hexagonal 
HoMnOa^ there is a transition between two symmetry 
nonequivalent antiferromagnetic phases. Wide domain 
walls where the spins rotate smoothly appear close to 
this transition, producing a pronounced magnetoelectric 
behavior. Interestingly, this behavior is not observed as- 
sociated with the abrupt domain walls where magnetiza- 
tion reversal is coupled to a change in the ferroelectric or- 
der parameter. In helicoidal magnetic manganites, with 
a perovskite crystal structure, the study of the dielectric 
dispersion and its behavior varying the temperature^S al- 
lowed to identify the motion of domain walls as the origin 
of the substantial enhancement of the dielectric constant 
in materials like DyMnOs and TbMnOs. 

Of particular motivational interest for the work de- 
scribed here is Ref. |29|, where manipulation and charac- 
terization of domain walls in ferroelectric and antiferro- 
magnetic BiFeOs was demonstrated. The authors of that 
effort could create ferroelectric domain walls and showed 
that domain walls with particular orientations exhibited 
larger conductance than others, within an otherwise good 
insulating material. Contrary to previous work, the au- 
thors focused on electronic transport in the direction par- 
allel to the domain walls, instead of transport across them 
as in other efforts in manganites »^^— A very recent work 
studied domain walls in a non-perovskite oxidei^ Do- 
main walls in hexagonal YMnOa have been shown to be 
more insulating than in the bulk, an effect actually oppo- 
site to the BiFeOs case. Here, motivated by those previ- 
ous efforts, we study the electronic properties of domain 
walls in undoped perovskite manganites. These materials 
are not only antiferromagnetic insulators, but also have 
a phase transition above room temperature, where an or- 
dered pattern of distorted MnOg octahedra sets in. We 
show that at finite temperature (T), structural domain 
walls display a conductance orders of magnitude larger 
that the conductance observed in bulk materials. 
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FIG. 1: Jahn- Teller domain wall at r=0, as obtained by the 
classical variables minimization process described in the text. 
The cluster size used is 4x28. (a) Variation of the order pa- 
rameter perpendicular to the domain wall, (black squares are 
results for A=1.25; red circles are for A=l). (b) Geometry for 
the calculation of the conductance, see text, (c) 2x 12 portion 
of the simulation cluster near a domain wall showing octahe- 
dra distortions in real space, for A=1.25. The oxygen displace- 
ments obtained by energy minimization have been rescaled to 
clearly show the two different lattice configurations left and 
right of the domain wall. 



In particular, we concentrate on the insulating, Jahn 
Teller (JT) ordered phase of undoped manganites with 
an A-type spin antiferromagnetic orderi^i^^ Our stan- 
dard variational and Monte Carlo calculations, using 
well-tested and reliable models, predict that: (i) domain 
wall electronic states exist, and are localized in the direc- 
tion perpendicular to the domain wall. They are anal- 
ogous to surface states and appear associated with the 
presence of an structural domain wall, (ii) These states 
form narrow bands within the energy gap of the bulk 
material. The periodicity of the structural distortions 
determine the position and width of theses bands. (Hi) 
At low temperature T there is a gap between the domain 
wall bands, and as T increases, conformational changes 
in the domain wall make the localized bands wider, and 
induce some spectral weight at the Fermi energy, leading 
to a finite conductance in the direction parallel to the 
domain walls. 



II. MODEL AND TECHNIQUES 
A. Model Hamiltonian 



The Hamiltonian used here is given by: 



H ^Hi 



DE 



Hr 



-ph 



(1) 



i^DE is the well-known two-orbital double exchange 
model Hamiltonian. This model describes the kinetic 
energy of the Cg electrons and their interaction with the 



magnetic background of the i2g core spins. ^"^ More pre- 
cisely, it is given by 



i/DE-- E A.-*7,7'4 



-/'^jn' 



(2) 



Here c| creates an electron at the Mn site i, in the 
eg orbital 7 (7 = 1,2 with 1 = \x'^ - y^) and 2 = 
|3^2 _ |.2^-j_ rp-j^g factor /^^ affecting the hopping am- 
plitudes in the limit of an infinite Hund's coupling de- 
pends on the Mn core spins orientation given by the an- 
gles 9 and ^ via the double-exchange mechanism /i,j = 
cos(6',/2) cos(6'j/2) + exp[i(V'i - Vi)] sin(6li/2) sin(6lj/2). The 
actual hopping amplitudes t take into account the dif- 
ferent overlaps between the orbitals along the directions 



n=x,y 



34 



1 1 — "^ 2 2 — V '^^1 ' 



-t and t\ 1=3^2 2~ 



-V^th 



--t. 



t is taken as the energy unit throughout this work, and 
its value depends on material details. Comparison be- 
tween theoretical calculations and experiments suggests 
that its value is of the order of half an electron volt J^ 
The phononic portion of the Hamiltonian reads: 

Ho-ph — ^2^ { — QliPi + Q2iTxi + QsiTzi) 

i 



where pi is the density operator at site i, and r^i and Tzi 
are the corresponding Pauli matrices in the eg subspace, 
that express the coupling of the lattice distortions, Q, to 
the electrons. In particular, Qu is the breathing mode 
of the MnOg octahedra around the z-th manganese ion, 
and Q2i, Qsi are the Jahn Teller modesi^ /3 = 2 is the 
stiffness of the breathing mode that effectively takes into 
account suppressed charge fiuctuations due to electron- 
electron interactions)^ 



B. Calculation Method 

Undoped manganites, such as LaMnOa, crystallize in 
an orthorhombic structure, with two nonequivalent man- 
ganese positions. Hamiltonian ([1]) has cubic symme- 
try, but in order to better compare with experiments 
we maintain the axis and notation corresponding to the 
orthorhombic unit cell, in particular we call b axis the 
direction perpendicular to the domain walls (see Fig[T]). 
We focus on the antiferromagnetic ^-type phasej-^? ex- 
perimentally found for La manganite and other trivalent 
cations of similar size.'^^ In this A-phase there are fer- 
romagnetic planes coupled antiferromagnetically in the 
c-axis direction. Double exchange effectively decouples 
electrons in the ab atomic planes, and 2D systems pro- 
vide a good description of the physics of these materials. 
In real compounds, certainly a coupling exists along the c 
axis that prevents 2D fluctuations to dominate the mag- 
netic and electronic properties of these compounds. 
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FIG. 2: Conductance as a function of temperature for dif- 
ferent system sizes. DW stands for conductance along do- 
main walls (see Fig. [TJ and Bulk, for conductance in the bulk 
limit. Details about the Monte Carlo procedure are provided 
in the text. Despite some small size effects, the conductance 
is clearly qualitative different when domain walls are present. 
For the bulk case, data for all sizes fall within the same curve 
in the scale of the figure. Note that the conductance is roughly 
independent of the dimensions along the b axis (perpendicular 
to the walls) , consistent with the hypothesis that the conduc- 
tance takes place along the domain walls. Also note that as- 
suming a hopping amplitude t of approximately 0.5 eV, then 
T = 0.06t is approximately room T. 



The Hamiltonian (IlJ depends on the octahedra and 
core spin configurations, Qi and Si. These have been 
approximated as classical degrees of freedom, as in most 
of the theoretical literature on manganites, and they de- 
termine the electronic properties of the system. At finite 
T, Monte Carlo simulations allow us to obtain the rel- 
evant equilibrium configurations and calculate thermal 
averages. Qi and Si are used as the Monte Carlo vari- 
ables and the resulting quadratic Hamiltonian for each 
configuration is exactly diagonalized using standard nu- 
merical routines)^ Further details about the Monte Carlo 
method and its application to manganites can be found 
in Refs. [33 andlsg. 

At low enough T, only one classical configuration be- 
comes relevant, and it can be determined by a suitable 
minimization algorithm. The results for T=0 presented 
here have been obtained with the Broyden method4S Al- 
though it still involves the diagonalizing of the fermionic 
sector for each step, this minimization process is faster 
than the standard Monte Carlo method and it can be 
applied to calculate results on larger clusters. Com- 
parison between results obtained with the two methods 
shows that size effects are small. This calculation method 
has been successfully applied to study uniform phases in 
manganites i^i^i^ Here, with the appropriate boundary 
conditions discussed below, we study the properties of 
structural domain walls in undoped Manganites. 



The ground state configuration for Hamiltonian ([T]) 
with one electron per site is a spin ferromagnetic plane 
with an ordered pattern of Q2 distortions characterized 
by a (tt, tt) wave vector in the 2D cubic notationi^i^ 
These orderings correspond in 3D to the experimentally 
observed A-type antiferromagnetic phase, and the (tt, tt, 
0) JT ordered phase. '^^ In the orthorhombic notation, the 
two octahedra around the two Mn ions in the unit cell 
have opposite Q2 distortions, one is elongated along the 
X axis, and the other along the y axis. The two possible 
ground state configurations intrinsic of a staggered or- 
der parameter (one reachable from the other by merely a 
global shift in one lattice spacing in any axes direction) 
are obtained by repeating the two possible units in the 
a and b directions. We label Q2 "^^ the order parame- 
ter, whose value is the magnitude of the distortions and 
its sign distinguishes between the two possible configura- 
tions. 

Monte Carlo or optimization calculations at low tem- 
peratures with standard periodic boundary conditions 
converge to one of the two possible ground state con- 
figurations. In order to study the domain walls that we 
wish to focus on in this effort different boundary condi- 
tions have to be imposed. In a MxN cluster, we fixed 
sites with coordinate 6=1 and b—N to the bulk equilib- 
rium values corresponding to Q2 "^^ >0, while sites at 
the center layer b=N/2 are constrained to have distor- 
tions with the same magnitude but corresponding to the 
other configuration Q2 "^^ <0. This way, two JT domain 
walls do appear in the system, and periodic boundary 
conditions can still be used for the fermionic sector. 

The high computational cost of the repeating diago- 
nalizing step of the Monte Carlo algorithm limits the 
sizes of the simulation cells. The largest cell used in this 
work is a=4, &=20, which contains 160 Mn ions (there 
are two Mns in the orthorhombic unit cell) . Thermaliza- 
tion is achieved by 2,000 Monte Carlo steps (as defined in 
Ref.[33) and measurement averages of 3,000 Monte Carlo 
configurations are used, taken every three steps to reduce 
self-correlations. This computational effort corresponds 
to roughly two days of running in a typical workstation 
for a 4x16 system and five days for a 4x20 system, for 
each fixed set of parameters in the problem. 

The conductance is calculated within the Landauer for- 
malism, as explained in Ref. |43J. To minimize size effects, 
we calculate the Green function of the system attached to 
leads that are identical copies of itself. Leads are enlarged 
by adding more copies of the system until the Green func- 
tion converges, providing an efficient method to calculate 
the intrinsic conductance)^ As illustrated in Fig. [1] the 
conductance of domain walls refers to conductance eval- 
uated in the direction parallel to the domain walls. 



III. RESULTS 

Before addressing the results regarding electronic 
transport, let us discuss the nature of the Jahn Teller 
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FIG. 3: Density of states a T—Q. States appear inside the 
bulk gap energy region when close to a domain wall (more 
specifically DW refers to the projected DOS at the line of 
manganese sites the closest to the domain wall in Fig. 1(c)). 
The states in this band are localized in the direction perpen- 
dicular to the domain wall (their weight decays exponentially 
in that direction) , but a small dispersion exists parallel to the 
domain wall. The band arising from the two eg orbitals is 
split due to the periodicity of the JT distortions in the direc- 
tion parallel to the domain wall. For A=1.25 the gap is of the 
order of t («0.5 eV), even in the presence of a domain wall. 
For the smoother domain wall corresponding to A=l the gap 
is very small and it cannot be resolved in this scale (inset). 



domain walls. Figure [Ija) presents the variation of the 
order parameter as a function of distance across the do- 
main wall. In magnetic domain walls, it is well known^ 
that the width varies as y^D/A, where A is the mag- 
netic anisotropy and D is the spin-stiffness. In the case 
of an isolated Jahn Teller center, the energy depends on 
{Q2+Q3) and the anisotropy is zero. The cooperative 
nature of Jahn Teller effect in manganites and the com- 
petition with the kinetic energy double-exchange term 
(even in the insulating state) gives rise to an effective 
anisotropy that depends on the value of A. These same 
effects stabilize the Q2 order in the ab plane.— Therefore, 
X/t controls the thickness of the domain wall. Figure [T] 
(a) shows how an abrupt domain wall is the minimum 
energy configuration for A=1.25, but for A=l, the do- 
main wall expands over several unit cells. The first A 
value provides a more accurate description of the real 
materials, as it can be deduced by the magnitude of the 
electronic gap discussed later. A real-space view of the 
octahedra distortions in the ab plane for the A=1.25 case 
is presented in Fig. [Ijc). 

The conductance vs. T plots contained in Fig. [2] con- 
stitute the main results of this manuscript. The con- 
ductance, G, refers to the conductance in the direction 
parallel to the domain wall, as shown in Fig. [Ijb) , and it 
is provided as a function of T. These results have been 
obtained by means of Monte Carlo simulations with dif- 
ferent system sizes, as indicated in Fig. [5] G is essentially 
zero (within our numerical precision) for bulk configura- 




FIG. 4: JT domain wall profile for r=0.01 and T=0.06 ob- 
tained via Monte Carlo simulations, at A=1.25 and for an 
a=3 6=20 system. The error bars show standard deviation 
due to thermal fluctuations. The temperature makes the do- 
main walls smoother, reducing the energy gap. 



tions, when no domain walls are present, for all the Ts 
considered in Fig. [21 However, when domain walls are 
present in the system, the conductance shows a clear up- 
turn when increasing the temperature. For A=l, where 
the domain wall is fairly wide, this upturn takes place 
at T sa0.02, while a larger T wO.OS is needed to induce 
a significant value for G in the A=1.25 case. There are 
some size effects in the results, but the increase with 
temperature of the conductance along domain walls, as 
compared to its negligible bulk value, is clearly present 
for all system sizes within our computational capabilities. 

Figure [3] shows the density of states (DOS) projected 
over different sites at zero T. It has been obtained by 
providing to the eigenenergies arising from the exact di- 
agonalization procedure a small lorentzian width (0.01 
t). For the bulk case (no domain walls) the DOS has 
a large gap due to Jahn Teller effects. This gap scales 
with A: it is roughly 2i for A=1.25 and half that value for 
A=l. Considering that t is in the range 0.5-1 eV, A=1.25 
is more consistent with the experimental estimations of 
the gap (1.6 eV in Ref. \M and 2.5 eV in Ref. \M)- 

When the density of states is projected over sites fur- 
ther than « 4 unit cells from the domain wall, we re- 
cover essentially the bulk density of states. However, a 
projection over Mn sites close to the domain wall is more 
interesting. Figure [3] shows that the DOS the closest to 
the domain wall has some features that do not appear in 
the bulk region. Two narrow peaks appear at energies 
located within the bulk energy gap. These two peaks are 
separated in energy creating a new, smaller, gap that de- 
pends strongly on A. It is sa i/10 for A=l and increases 
to ~ i for A=1.25. The integrated weight of these peaks 
decays exponentially in the b direction, perpendicular to 
the domain wall, so these states are localized at the do- 
main walls. A least-squares fitting shows that the decay 
length is of the order of the nearest neighbors distance. 
The fermionic wave functions coming out of the simu- 



lations, via the exact diagonalization procedure, have a 
weh defined momentum in the a direction. Therefore, 
we can state that the new peaks correspond to a band 
locahzed in the b direction and with a small dispersion 
(« i/20) parallel to the domain wall. Note that the arti- 
ficial broadening 0.01 1 is not sufficient to understand the 
width of the in-gap peaks; in fact, the actual dispersion 
in the a direction produces their intrinsic width. 

These bands are obvious candidates to explain the con- 
ductance results in Fig. [51 Although the energy gap due 
to the splitting of the localized band is rather small for 
A=l, it is too large to explain the nonzero conductance at 
T=0.05 observed for the more realistic case A=1.25. As 
the electronic states associated with the domain wall dis- 
perse in the a direction, some splitting due to the periodic 
Jahn Teller distortions in that direction can be expected. 
This explains why the splitting decreases with A. Note 
that the conductance calculation has enough precision to 
determine that the small gap in the A=l case, difficult 
to appreciate in Fig. [3l still suppresses transport at the 
lowest temperatures. 

In order to understand the onset of electronic trans- 
port with temperature, we examine the structure of the 
domain walls as T increases (Fig. U]), for A=1.25. We 
have observed that the order parameter changes much 
more abruptly at T=0.01 as compared to T=0.06. The 
points in Fig. [3] correspond to average values over 1,500 
Monte Carlo steps, with the error bars showing the stan- 
dard effect of thermal fiuctuations. The domain wall is 
clearly wider for r=0.06, with a width that effectively 
corresponds to a smaller A if the calculation were done 
at very low T (as discussed below). Entropy reduces 
the effective anisotropy (or lowers the order parameter 
stiffness). A visual inspection of Monte Carlo snapshots, 
namely equilibrated configurations of classical variables, 
shows that the widening of the domain wall is due to 
changes in the intrinsic width, and the contribution from 
thermal vibrations of the center of the wall is small. 

The widening of the domain wall induces changes in the 
density of states. In Fig. [5] we show the DOS projected 
on sites at different distances from the domain wall, ther- 
mally averaged at T=0.06. Particle-hole symmetry was 
assumed, and the spectrum symmetrized with respect to 
Ey = in order to reduce the noise (Monte Carlo calcu- 
lations are more noisy than optimization results due to 
thermal and finite size effects). However, it can be clearly 
observed how the band associated with the domain wall 
is now wider in energy, and has some spectral weight at 
the Fermi energy, that is responsible for the finite con- 
ductance. Notice that this band is still localized in the 
h direction and by moving 4 sites away from the domain 
wall, the bulk density of states is recovered. 

An intuitive way to understand the results presented 
above is to consider that the coupling A is now effectively 
smaller near the domain wall. But note that temperature 
must also be incorporated. At T=0, the domain wall is 
abrupt for the realistic coupling A=1.25, and therefore 
this coupling produces distortions at the domain wall as 
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FIG. 5: Thermal averaged electronic density of states for 
r=0.06, A = 1.25, projected over different atomic planes par- 
allel to the domain wall (in practice these are atomic chains 
since we use a 2D geometry). The zero of reference for each 
case is provided by the result at the far left or right in the fig- 
ure. At finite temperature the domain walls are not as abrupt 
as at r=0, the bands associated with it (see Fig. [3]) are wider 
and some spectral weight is shifted to the Fermi energy _Ef. 
The eigenenergies arising from the Monet Carlo simulation 
have been given a small width (O.Olt), and the spectrum has 
been symmetrized with respect to E-p to reduce noise. 



large as in the rest of the system. The bands appear 
associated with the domain wall due to the change in 
periodicity. For smaller A's, domain walls are smoother. 
Wider domain walls have a smaller gap, relative to the 
bulk gap. For A=1.25, the gap between the domain walls 
gaps is half the bulk value, while is an order of magnitude 
smaller than the bulk gap for A=l. Thus, T must also be 
included to justify the effective reduction of the value of 
the electron lattice coupling. From the structural point 
of view, this effect results in a smaller anisotropy and an 
increase in the domain wall width. As regarding the elec- 
tronic structure, it reduces the gap between the localized 
bands, eventually closing it. 



IV. CONCLUSIONS 

Stimulated by the recent experimental characterization 
and manipulation of structural domain walls in a per- 
ovskite transition metal oxide, ^^ we have undertaken the 
theoretical study of domain walls in undoped manganites. 
These materials display an ordered phase with an alter- 
nating pattern of distorted MnOg octahedra, and mobile 
eg electrons are known to couple strongly to these dis- 
tortions. The equilibrium configurations and electronic 
properties of structural domain walls have been exam- 
ined, and they remarkably differ form the properties of 
the bulk. Since it is an important characteristic for ap- 
plications and a property directly measurable by exper- 
iments, here special attention has been paid to the elec- 
tronic conductance of the system. 

Electronic transport along structural domain walls is 



dominated by two narrow electronic bands associated 
with them. These domain walls electronic states have 
energies within the bulk energy gap region. At low T's, 
domain walls are abrupt, and the narrow bands are sepa- 
rated in energy. The energy gap between them, although 
smaller than the bulk energy gap, is still large enough to 
make the conductance negligible at low T. However, for 
T's of the order of room temperature, thermal fluctua- 
tions favor wider domain walls. This makes the disper- 
sion of the domain walls electronic states larger, reduc- 
ing the energy gap, and inducing some spectral weight at 
the Fermi energy. Correspondingly, conductance along 
domain walls is several orders of magnitude larger for 
the temperature range studied than in the bulk, while it 
is still zero within our numerical precision for the bulk 
phase. Our results pose the question of whether the con- 
ductance enhancement by domain walls electronic states 
induced by changes in the periodicity of the system might 



be a more general phenomenon, and suggest the possi- 
bility that they might play a role in recent experiments 
showing conducting domain walls in a multiferroic per- 
ovskite oxide. ^^ These issues certainly deserve further 
theoretical and experimental studies. 



Acknowledgments 

The authors acknowledge useful conversations with 
Jan Seidel. This work was supported by the NSF grant 
DMR-0706020 and the Division of Materials Science and 
Engineering, Office of Basice Energy Sciences, U.S. De- 
partment of Energy. R.Y. acknowledges support from 
NSF Grant No. DMR-0706625, the Robert A. Welch 
Foundation Grant No. C-1411, and the W. M. Keck 
Foundation. 



^ E. Dagotto, Science 309, 257 (2005). 

^ S. Cheong and M. Mostovoy, Nature Materials 6, 13 

(2007). 
^ W. Eerenstein, N. Mathur, and J. Scott, Nature 442, 759 

(2006). 
'^ C. Kittel, Rev. Mod. Phys. 21, 541 (1949). 
^ E. Lifshitz, L. Landau, and L. Pitaevskii, Electrodynamics 

of continuous media (Pergamon Press, 1984). 
® C. Marrows, Adv. Phys. 54, 585 (2005). 
'' G. Tatara and H. Kohno, Phys. Rev. Lett. 92, 086601 

(2004). 
® M. Klaui, P.-O. Jubert, R. Allenspach, A. Bischof, J. A. C. 

Bland, G. Faini, U. Riidiger, C. A. F. Vaz, L. Vila, and 

C. Vouille, Phys. Rev. Lett. 95, 026601 (2005). 

^ G. Tatara and H. Fukuyama, Phys. Rev. Lett. 78, 3773 

(1997). 
^° P. Schiffer, A. P. Ramirez, W. Bao, and S.-W. Cheong, 

Phys. Rev. Lett. 75, 3336 (1995). 
" D. L Golosov, Phys. Rev. B 67, 064404 (2003). 
^^ Q. Li, Y. Hu, and H. Wang, Journal of Applied Physics 

89, 6952 (2001). 
" P. M. Levy and S. Zhang, Phys. Rev. Lett. 79, 5110 (1997). 
^^ N. Mathur, P. Littlewood, N. Todd, S. Isaac, B. Teo, 

D. Kang, E. Tarte, Z. Barber, J. Evetts, and M. Blamire, 
Journal of Applied Physics 86, 6287 (1999). 

^5 Y. Suzuki, Y. Wu, J. Yu, U. Ruediger, A. Kent, T. Nath, 
and C. Eom, Journal of Applied Physics 87, 6746 (2000). 

^® M. Yamanaka and N. Nagaosa, Physica B: Condensed Mat- 
ter 237, 28 (1997). 

^'' T. Becker, C. Streng, Y. Luo, V. Moshnyaga, B. Dam- 
aschke, N. Shannon, and K. Samwer, Phys. Rev. Lett. 89, 
237203 (2002). 

^* G. Singh-Bhalla, S. Selcuk, T. Dhakal, A. Biswas, and 
A. F. Hebard, Phys. Rev. Lett. 102, 077205 (2009). 

^^ T. Arnal, A. V. Khvalkovskii, M. Bibes, B. Mercey, 
P. Lecoeur, and A.-M. Haghiri-Gosnet, Phys. Rev. B 75, 
220409(R) (2007). 

^° V. L. Ginzburg and L. D. Landau, h. Eksp. Teor. Fiz. 20, 
1064 (1950). 

^^ S. K. Streiffer, J. A. Eastman, D. D. Fong, C. Thompson, 



A. Munkholm, M. V. Ramana Murty, O. Auciello, G. R. 

Bai, and G. B. Stephenson, Phys. Rev. Lett. 89, 067601 

(2002). 
^^ P. Aguado-Puente and J. Junquera, Phys. Rev. Lett. 100, 

177601 (2008). 
^■^ J. Pi^ivratska and V. Janovec, Ferroelectrics 204, 321 

(1997). 
^^ A. Aird and E. K. H. Salje, Journal of Physics: Condensed 

Matter 10, L377 (1998). 
^^ M. Bartels, V. Hagen, M. Burianek, M. Getzlaff, U. Bis- 

mayer, and R. Wiesendanger, Journal of Physics: Con- 
densed Matter 15, 957 (2003). 
^^ M. Fiebig, T. Lottermoser, D. Frqohlich, A. Goltsev, and 

R. Pisarev, Nature 419, 818 (2002). 
^'^ T. Lottermoser and M. Fiebig, Phys. Rev. B 70, 220407(R) 

(2004). 
^** F. Kagawa, M. Mochizuki, Y. Onose, H. Murakawa, 

Y. Kaneko, N. Furukawa, and Y. Tokura, Phys. Rev. Lett. 

102, 057604 (2009). 
^^ J. Seidel, L. Martin, Q. He, Q. Zhan, Y. Chu, A. Rother, 

M. Hawkridge, P. Maksymovych, P. Yu, M. Gajek, et al.. 

Nature Materials 8, 229 (2009). 
^° T. Choi, Y. Horibe, H. Yi, Y. Choi, W. Wu, and S. Cheong, 

Nature Materials 9, 253 (2010). 
^^ T. Hotta, M. Moraghebi, A. Feiguin, A. Moreo, S. Yunoki, 

and E. Dagotto, Phys. Rev. Lett. 90, 247203 (2003). 
^^ J. Rodn'guez-Carvajal, M. Hennion, F. Moussa, A. H. 

Moudden, L. Pinsard, and A. Revcolevschi, Phys. Rev. B 

57, R3189 (1998). 
^^ E. Dagotto, Nanoscale Phase Separation and Colossal 

Magnetoresistance (Springer- Verlag, Berlin, 2002). 
^'^ J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 
^^ J. Salafranca, G. Alvarez, and E. Dagotto, Phys. Rev. B 

80, 155133 (2009). 
^^ T. Hotta, S. Yunoki, M. Mayr, and E. Dagotto, Phys. Rev. 

B 60, R15009 (1999). 
^■^ E. WoUan and W. Koehler, Phys. Rev. 100, 545 (1955). 
^* E. Anderson, Z. Bai, C. Bischof, S. Blackford, J. Dem- 

mel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hanimar- 

ling, A. McKenney, et al., LAPACK Users' Guide (Society 



for Industrial and Applied Mathematics, Philadelphia, PA, 
1999), 3rd ed. 

^^ J. Salafranca and J. A. Verges, Phys. Rev. B 74, 184428 
(2006). 

■"^ W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical recipes in FORTRAN, The art of sci- 
entific computing (1992). 

■'^ C. §en, G. Alvarez, and E. Dagotto, Phys. Rev. B 70, 
064428 (2004). 

■'^ J. Salafranca and L. Brey, Phys. Rev. B 73, 024422 (2006). 



^^ J. Verges, Comput. Phys. Commun. 118, 71 (1999). 

'*'' Z. Popovic and S. Satpathy, Phys. Rev. Lett. 84, 1603 

(2000). 
''^ K. Murakami, T. Yamauchi, A. Nakamura, Y. Moritomo, 

H. Tanaka, and T. Kawai, Phys. Rev. B 73, 180403(R) 

(2006). 
■'^ T. Inami, T. Fukuda, J. Mizuki, S. Ishihara, H. Kondo, 

H. Nakao, T. Matsumura, K. Hirota, Y. Murakami, 

S. Maekawa, et al., Phys. Rev. B 67, 045108 (2003). 



